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We investigate an analytically tractable toy model for thermally induced polymorphic dynamics of 
cooperatively rearranging biofilaments - like microtubules. The proposed 4 -block model, which can 
be seen as a coarse-grained approximation of the full polymorphic tube model, permits a complete 
analytical treatment of all thermodynamic properties including correlation functions and angular 
fourier mode distributions. Due to its mathematical tractability the model straightforwardly leads 
to some physical insights in recently discussed phenomena like the "length dependent persistence 
length". We show that a polymorphic filament can disguise itself as a classical wormlike chain on 
small and on large scales and yet display distinct anomalous tell-tale features indicating an inner 
switching dynamics on intermediate length scales. 

PACS numbers: 87.16.Ka, 82.35.Pq, 87.15.-v 



INTRODUCTION 

Biological filaments of the cytoskeleton are large 
macromolecules formed by a hierarchical self-assembly of 
smaller yet often highly complex protein subunits. The 
monomer complexity can allow the subunits to undergo 
rearrangements between several conformational states. 
Once inserted into a macromolecular lattice these indi- 
vidual subunits can start to interact giving rise to co- 
operative phenomena which can affect equilibrium and 
dynamical properties of the whole assembly in unex- 
pected manners. It is the goal of this paper to explore 
in detail the statistical mechanics of this type of co- 
operatively switching supramolecular assemblies, whose 
paradigm example might be found in the stiffest cy- 
toskeletal protein filaments of eukaryotic cells, the mi- 
crotubules [l[ . Microtubules are hollow nanotubes whose 
walls are formed by lateral self-association of parallel 
protofilaments that themselves are built by a head-to-tail 
polymerization of a/3-tubulin heterodimer protein sub- 
units. This very complex architecture confers to micro- 
tubules their high stiffness as well as a number of unique 
static and dynamic properties. In Q, guided by exper- 
imental findings, we built the case for a novel model of 
microtubules, based on two hypotheses: the elementary 
tubulin dimer units can fiuctuate between a curved and 
a straight configuration and can interact cooperatively 
along a protofilament axis. This implies that the ground 
state of microtubules is not, as usually accepted a straight 
Euler beam, but instead a fluctuating several micron 
sized cooperative super-helix. The resulting polymor- 
phic dynamics of the microtubule lattice seems to quan- 
titatively explain several experimental puzzles including 
anomalous scaling of dynamic fluctuations of grafted mi- 
crotubules 0, their apparent length-stiffness relation Q 
and their remarkably curved-helical appearance in gen- 
eral These results rely on phenomenological mod- 
elling where the cylindrical symmetry of the microtubule 



lattice is approximated by a continuous symmetry [2j. 
This approximation seems reasonable as the number of 
protofllaments is typically large, ranging from 9 to 17 for 
taxol stabilized in-vitro microtubules with a predominant 
14 protofllament structure [1,0. In vivo, microtubules 
most commonly appear with 13 protofilaments 's'l, al- 
though many exceptions exist depending of the cell type. 
The approach developed in Q revealed the existence of 
a unique and unusual zero mode dynamics which has 
strong consequences : microtubules are permanently co- 
herently reshaping -i.e. changing their reference ground 
state configuration- by thermal fiuctuations. 

Going beyond the continuous phenomenological ap- 
proach of [2], we explore here other important aspects 
of the polymorphic microtubules theory by considering 
a simplified minimal solvable model. We will adopt a 
coarse-grained approach where the microtubule lattice is 
considered as made up of only 4 blocks of protofilaments 
that can fiuctuate between a curved and a straight con- 
figuration. Although the coarse-grained model to some 
extent loses the (quasi)continuous zero mode dynamics, 
it captures a number of important features of the full 
model [il and exact analytic computations of relevant 
observables becomes accessible in this simpler case. In 
particular the tubulin-tubulin state correlation function, 
the persistence length, the thermodynamic stability of 
conformational states as well as the tangent angle spec- 
trum can be computed rather elegantly. This simplified 
approach should give more detailed analytical insights 
into microtubule's static properties and make connection 
with the previous phenomenological model's parameters 
introduced in Q. While it is not our prime goal here to 
compare the developed toy model with experiments we 
hope that some of the derived results (like the tangent an- 
gle power spectrum) can become useful guides for future 
experimental quantiflcation of polymorphic fllament fluc- 
tuations. For a more detailed experiment-theory com- 
parison and deeper motivation of the model we refer the 
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reader to the articles 0. 



THE 4-BLOCK POLYMORPHIC TUBE MODEL 

The microtubule lattice (see Fig. 1) is modelled as 
a continuum tube material made of a variable number 
of protofilaments with Ri « 7.5nm, Ro ~ 11.5nm - 
the inner and outer microtubule radii respectively. The 
protofilaments are twisted around the microtubule's lon- 
gitudinal axis with the corresponding internal twist qq 
- or equivalently the pitch A = inq^^ - being a lat- 
tice type dependent constant that takes typical discrete 
values A — -1-3.4/im, +25fim, ~6/im for 12, 13 and 14 
protofilament microtubules respectively d, Q [1] • 

To describe the tube's geometry we introduce two ref- 
erence frames (cf. Fig. lA). One is the material frame 
with base vectors {ui,U2,uz) attached to each micro- 
tubule cross-section. The other is an external fixed lab- 
oratory frame with base vectors {ux, Uy,Uz) with respect 
to which all deformations are measured. Putting the fil- 
ament along the Uz axis direction and considering small 
angular deflections we have Uz ~ u^- In this case the two 
frames are simply related by a rotation transformation 
i?(s) given by internal microtubule lattice twist go: such 
that {ux,Uy) = R{s){ui,U2) with 
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and s G [0, L] the longitudinal position variable along 
the microtubule centerline. The arbitrary rotation angle 
if corresponds to the angular deviation between the two 
frames at s — 0. As each protofilament is built by a self 
association of N discrete GDP-tubulin dimers of length 
b ss 8nm it is natural to introduce a discrete variable 
i such that all microtubule's cross-section positions are 
written as s = ib with i = 1..N. 

The coarse-grained block approximation consists of 
grouping neighboring protofilaments into n — 4 effec- 
tive blocks (cf. Fig. IB). The fiuctuations of the block- 
dimers between 2 states - a straight and a curved state 
with intrinsic curvature kfp and with an energy differ- 
ence AG > favoring the curved state are modelled by a 
two state variable trf = 0, 1 where i = 1..N is the longi- 
tudinal position and k — 1..4 the block's index (cf. Fig. 

IC) . 

To complete the description of the polymorphic tube 
model two order parameters at each microtubule cross- 
section can be introduced [2)]. The first is the vectorial 
polymorphic order parameter pi = wi(cr| ^crf)+ 'U2{o'f — 
af), a 2D vector attached to each cross-section (cf Fig. 

ID) describing the asymmetry of the curved state distri- 
bution - a kind of "conformational polarization vector" 
of the block states. For instance the "all-straight" or 
"all-curved" protofilament state correspond both to the 




FIG. 1: Elements of the "polymorphic four block tube 
model", (a) Tube contour with the internal and external 
frames, (b) A short tube section with intrinsic twist go with 
four blocks at each section, (c) The blocks can fluctuate co- 
operatively between two discrete states. The outwards curved 
a = 1 (arrow) state is energetically preferred over the straight 
state a — (circle) with an energy gain E = —AG. The junc- 
tion between straight and curved states along the longitudinal 
block axis are penalized by a coupling constant E = +J. (d) 
Competition between block switching energy and elastic lat- 
tice strain energy leads to spontaneous symmetry breaking: 
the tube bends to a randomly chosen direction and assumes 
a non-zero vectorial polymorphic order parameter p ("confor- 
mational polarization") (e) The four block model can be de- 
coupled into a pair of two-block models which can be mapped 
onto a ladder type Ising model with coupling constants — _B 
and +J and a bias field 7. 



same value p = 0, cf. Fig. ID (as the curved state dis- 
tribution is isotropic in both special cases). A second 
(scalar) quantity rrii = X)fc=i '^i counts the total number 
of blocks in the curved state at each cross-section i. The 
tubulin cooperativity along each protofilament (block) 
axis is modelled by an Ising type nearest-neighbor co- 
operative interaction with an interaction energy J > 
favoring longitudinal nearest neighbors to be in the same 
state. The interaction energy at cross-section i reads: 



Winter A 



(2) 
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The total elastic + polymorphic energy can then be writ- 
ten as 01 : 



EN 
(ei + einter,i) ■ 
z— 1 



(3) 



where is the combined elastic energy and the energy 
resulting from the switching of tubulin dimers at i-th 
cross-section: 
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with the elastic bending modulus B = ^ (i?^ — and 
Y the Young modulus. The effective lattice curvature 
Ki — ^"^^^^^2^ i^FP results directly from the preferred 

curvature of the individual protofilament upp Here 
K is the microtubule centerline curvature vector and Kpoi 
the polymorphic curvature vector which in the external 
coordinate frame {ux,Uy) is 
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In Eq. |4] we introduced an important dimensionless pa- 
rameter 7 = -^^^ — fp^A which measures the competi- 
tion between block switching and elastic energy and ulti- 
mately determines the microtubule shape. It effectively 
acts as an external field that biases the curved lattice 
state for 7 < and favors the straight state for 7 > 
(cf. below). For small deflections Or^/y << 1 around the 
z-axis, the unit vector tangent to the microtubule's cen- 
terline is approximately given by t « {9xT0y,l) in the 
laboratory frame {ux,Uy,Uz) where 9 = {0x,0y) are the 
centerline deflection angles in x/y direction. The global 
centerline curvature vector k = dt/ds can then be ap- 
proximated as w (O'xtO' ,0) . Writing the total cur- 

vature as 1^ = it pd + 6' ei with 6' the purely elastic 



contribution we readily see from Eqs. I2l3l that the parti- 
tion function decomposes into a product of two indepen- 
dent elastic and polymorphic contributions Z = Z^iZpoi, 

with Zei = / D6'exp (^-^-f Jo Ofids^ (and h = B/kpT 
the bending persistence length with kgT the thermal en- 
ergy). Therefore our main goal reduces to the computa- 
tion of Zpoi — J2a^^Pi~k^^cr) with Eg. the polymor- 
phic contribution of Emt- In the general case (with large 
number of blocks) this appears to be a formidable task, 
that however can be exactly performed in the case of 4 
blocks as we will see in the following. 

Reduction to a 2 x 2-block model 

We first note that, at each cross-section p| = 

{af — afY + i'^i ~ (^lY decomposes into a sum over mu- 
tually facing blocks i.e. only blocks 1-2 and 3-4 couple 
directly. Consequently the partition function can be fac- 
torized into a product of two simpler ones Zpd — ^12^34 



with Zki the partition function of the Ising variables 
and cr' of the blocks that face each other, cf. Fig. ID. 
Therefore from now on we may consider the 2 block 
model, say the block pair 1 — 2. The polymorphic or- 
der parameters are in this case mi — a} -\- and 
J^i = (a} — af) ui. By introducing new more convenient 
Ising variables — ±1 defined as S'f — 2crf — 1 the 
energy Eq. [2] takes the familiar form of a ladder type 
Ising model (cf. Fig. IE). In terms of these variables, 
Z12 = Z^Z with Zo = exp(^BKfL (717 - 1) /4) and 



E exp(-i?) 

{Si,S2} = -l,l 



(6) 



where from Eqs. [5]|3]we obtain 
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with B 



J 



iksT' " ~ ^^'^ ^ ~ fT-^- With the two 
spins - on either side of the ladder (Fig. IE)- being inter- 
changeable we have {S}) = (Sf) (and thus {a}) = (erf)) 
for all i as the energy is translationally invariant. From 
Eq. [5l we readily see that {Kpoi,i) = 0. And yet this 
does not mean that the microtubule is on average in 
a straight configuration. Just on the contrary, the mi- 
crotubule can form a three dimensional super helix Q 
that is "coherently" reshaping as the curvature alter- 
nately exchanges from one block to the other with a non- 
vanishing value of (^Kpoi^- Therefore the order parame- 
ter that characterizes the typical curvature of the lattice 
is the mean "polarization" P = ^/Tf^^^~Wf) of the 



and 7 



curved states (as ^n-po/ 
quantity is M = ;^ Ei (c^f ) 



The other important 

(cF^) - the mean concen- 
tration of curved states, or in the Ising terminology the 
mean "magnetization" (up to a trivial additive constant) 
. Within the same terminology the parameter 7 takes 
the role of a "magnetic field" that, according to its sign, 
favors one or the other possible spin orientation (curved 
or straight block) -but the same orientation for blocks on 
both sides of the ladder in Fig. IE. The "ferromagnetic 
coupling constant " J promotes the longitudinal parallel 
alignments of spins along a given block axis. The "anti- 
ferromagnetic coupling" B on the other hand favors S} 
and to be antiparallel - a tendency which competes 
with 7 and the lattice is frustrated. For large I7I >> B, 
we therefore expect the alignment tendency to win on 

both sides of the ladder so that l^S]^'^^ ~ ±1, (+1 if 

7 < and -1 if 7 > 0) and thus M « 1 or 0. In this sit- 
uation the microtubule is in a straight state with P « 0, 
which is either completely unstressed with M w (for 
positive 7) or maximally prestrained state with A/ ss 1 
(for negative 7). 

For 7 = 0, the lattice is not frustrated and on aver- 
age when on one ladder-side Sj = 1 then Sf = —1 on 
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the other (and vice versa). Consequently (Sf ) = and 
M = (erf) = 1/2. In this case blocks 1 and 2 at the 
cross-section i fluctuate alternately between the curved 
and straight state and the tube bends alternately in the 
directions Ml and —ui. The cooperative ferromagnetic in- 
teraction J implies a correlation between the spins along 
the contour and the formation of domains with size of the 
order of the spin-spin correlation length ^ computed be- 
low. It is these whole domains that alternately bend the 
tube in the Ui and —Ui direction that lead to {'1^ i) — 
but P 0. In fact P will take its largest value P w 1 in 
this non frustrated case. One remarks that for ^ much 
smaller that the internal twist wavelength A the typical 
domain looks like a circular arc, whereas for ^ 3> A, a 
typical coherent domain is a super-helix in the 3 dimen- 
sional space with pitch A as the direction of bending ±ui 
is slowly rotating in the external frame. For L >> ^, the 
microtubule is made up of a juxtaposition of uncorrelated 
fluctuating helices or of uncorrelated circular arcs that 
bend independently in the zLui direction. Consequently 
we expect a long distance behavior similar to an usual 
worm like chain. This is true despite the fact that elastic 
contributions where neglected in this qualitative discus- 
sion. Indeed there are no elastic fluctuations here (as 
for usual semiflexible filaments) : the polymorphic transi- 
tion from curved to straight states of short uncorrelated 
segments only mimic elastic fluctuations. For 7 small 
but non zero the lattice is slightly frustrated and thus 
less curved (smaller value of P) displaying nevertheless a 
qualitatively similar physical behavior. 

2-BLOCK MODEL THERMODYNAMICS 

The partition function Eq. [51 can be exactly computed 
via the transfer matrix method. For this purpose, we im- 
pose periodic boundary conditions on both sides of the 
ladder 5'^^ = S^^^. This permits us to write the parti- 
tion function as Z =Tt{T^) with the symmetric (4 x 4) 
matrix T defined by 

g - 1 ( + S + S 1+ 1 + S,"+ - 1 ( Si Si + S ■ + 1 1 ) + JS 1 S 1 + 15== S 1 

with explicit elements given by 

Denoting U the matrix diagonalizing T, such that 
UTU^^ = A is a diagonal matrix with eigenvalues A^ 
with r = 1,4 , the partition function can be written as 
Z = Tr{A^) = Er=i,4^i^- Denoting Ai the largest 



eigenvalue the free energy per lattice site f = F/N = 
— ksT/NlTiZ reduces in the thermodynamic limit of 
— ?> 00 to the expression / = — fc^rinAi. From the 
free energy, all other thermodynamic quantities can be 
derived. In particular the curved state density M = 

5 — I; ^ In Ai and the polarization P — ^\ -I- In Ai. 

The explicit expression for Ai is rather cumbersome 
and we omit it here. Instead we provide the plot of 
M and P in terms of 7 for different values of J with 
B — \ (in 0, J5 was indeed found to be of order unity) 
in Fig. 2. We first remark that P(— 7) = ^(7), so the 
microtubule's mean curvature is symmetric with respect 
to the sign of 7. Further we observe that there always ex- 
ists a range of 7 where P is close to unity which is quite 
spread for small coupling J <^ UbT and that for I7I ^ B, 
the polarization P is larger for smaller values of J. Al- 
though in this regime the total energy density is minimal 
for straight unstressed or prestrained states, these states 
have a small entropy. In comparison, curved tube states 
have a higher energy density in this regime but also a 
higher entropy, as blocks of small size ^ fluctuate inde- 
pendently. Consequently these curved -helical- states can 
have a smaller free energy. With growing coupling J, the 
size of coherent blocks increases (as ^ grows with J) and 
the energy contribution becomes dominant over entropy 
so P gets lower with growing J. From Fig. 2, we observe 
that for J = l.SfcsT, the entropy is already negligible as 
for I7I > B the straight states with P w are selected. 
On the contrary for I7I < B, the microtubule adopts a 
curved conformation with a quasi constant value P « 1- 
corresponding to maximum lattice curvature which is in- 
dependent of 7. This plateau region shows that in order 
to observe a curved (or helical) conformation microtubule 
does not require & very precise fine tuning of 7 - as long it 
is smaller than B. This approximate 7 independence al- 
lows us to limit ourselves to the analytically most elegant 
case 7 = 0. 



Simple case 7 = 

In this case all quantities can be computed in compact 
form. In particular the largest eigenvalue of the transfer 
matrix simplifies to Ai = 2 cosh(P) cosh(2 J) -I- D with 
D = (3 + cosh(2S) + 2cosh(4J)sinh^(S))i/2. As already 
mentioned the curved state density is constant M = 1/2 
(independent of B and J), whereas the total "polariza- 
tion" is P2 = 1 - 2 (alflf ) = i + ^i"h(5)^osh(2j) ^^.^j^ 

tends to 1 for J > 1. This means that {(j}<jf) which 
can be understood as both states (on opposing ladder 
sides) at the same site i fluctuate alternately between 
and 1. The coupling J which is a measure of the cooper- 
ativity determines over which distance the curved state 
is coherent. The cooperativity can be measured by the 
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FIG. 2: The curved state density M = ^ Y.'^^i {"i) ver- 
sus the state bias parameter 7 for different values of coupling 
constant J and fixed stiffness B = 1. Note the broad 7 inde- 
pendent plateau region for J = 1.5 where M = f/2. For large 
values of I7I , M is either or 1 corresponding to straight 
unstrained or straight strained lattice states respectively. 




the matrix U diagonalizing T, i.e. U ^TU = A and the 
following transformed matrices V'^ = U^^V^U , we have 
(S-fSj.) = iTr(A'-iy'=AJ-'y'A^-J+i). In the limit 
N ^ CO, this expression can be evaluated : 



A4 

Ai 



(8) 



^31^3 

where the eigenvalues are ordered Ai > .. > A4. From Eq. 
[8] all g^j can be deduced. It turns out that the most in- 
teresting quantity connected to the spatial fluctuations of 
the microtubule -its persistence length (see next section) 
- is connected to /^g,j = g}} - glf = (c^^^crj) - crj) . 
An explicit computation of the matrix elements of 
shows that V^^ = V^i = and V^iV^^ = Vi-^V^^ as well 
as V^^Vl^ = -Vf2^2\ = P'^- This leads to 



where the correlation length is ^ = 



(9) 



with 



A2 = 2exp(B)sinh(2J). With the result for 2 blocks 
(ladder model, Fig IE) we are now able to compute the 
persistence length for the 4-block polymorphic tube (Fig 
IB). 



ANGULAR PERSISTENCE LENGTH 




FIG. 3: Plot of the mean polymorphic polarization P of the 
reduced two block model versus the conformational bias field 
7 for different values of coupling constant J and effective tube 
stiffness B = 1. For J > 1.5, the tube's curvature is at maxi- 
mum for I7I < 1, as P « 1 in this region. For larger I7I , the 
tube is in a straight state and P goes to zero. 



spin-spin correlation function (7*!^? = (cr^^irj) — (trf ) (fij) 



(cr,fcrj) - 1/4 with k,l 
j. By symmetry gj^ = ^..^ y.^ - 
the two matrices and V'^ defined as 



1 , 2 at two different sites i and 
and gjf — gfl- Introducing 



/ 1 \ 
010 
00-10 
\ -1 / 



and = 



/ 1 \ 
0-100 
10 
\ -1 / 



the correlation functions of the variable 
with k ~ 1,2 can be written as (S^Sj) — 

^rj.^(^rj.,.-lykj.j-,^lylj.N-j+Z+r^ fOr i < j. With 



The persistence length Ip is the length scale char- 
acterizing the filament's resistance to thermally in- 
duced bending moments. For usual biofilaments like 
DNA or actin filaments Ip = is a material con- 
stant equal to the bending persistence length Ib = 
B/ksT. However for microtubules Ip is known to be 
length dependent Here we compute the angu- 

lar persistence length for the 4-block model from the 
usual definition in terms of the tangent-tangent cor- 
relation fmiction (l^ {s') ■ ~t (s)'^ = exp(— |s' — s| /Zp) 

where ~t{s) is the unit-tangent vector at the posi- 
tion s of the microtubule centerline. From t (s) = 
{0xis),8y{s), yjl — 8x{s) — ^'^(s)) in the external frame, 

we deduce to quadratic order in Q that {&') ■ ~t{s)'^ w 

1 — ^V{s' — s) with the angular variance V{s' ~ s) — 
{{OAs') ~ d-As)f) + {{0y{s') - Oy{s)f) . Therefore one 
can write ^~i^(s') • ~^(s)^ exp(— |s' — s| /Zp) with a 



persistence length lp{s' — s) = 2-^ 



.TTT-, — -\ that will be man- 

V(s'-s) 

ifestly distance dependent in our case. Now by writ- 
ing 6* = 6* pol 



6 el and from the independence of 
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polymorphic and elastic contributions (^9 poi ■ 6 

6 poi^ ■ (^0 ei^ =0, we deduce that the total persis- 
tence length can be decomposed as = (7^ + r~[)^^ 
with a polymorphic persistence length 



lpoi{s' - s) = 2 



Vpoi{s' - s) 



(10) 



At this level, for the geometric description of the mi- 
crotubule shape in space it is more convenient to use a 
continuum description and replace the discrete index i 
by the continuous variable s, so that 9 poi^i — >■ 9 poii-s) 
which is obtained from the integration over s of the cur- 
vature Eq. [5] 

^ f 
9 poi{s) ^ 9poi{0)+ / K.poi{s)ds (11) 
Jo 

Now we can compute the polymorphic vari- 
ance from the correlation function Gpoi.x{s',s) ~ 
{{OpoiAs') - OMiOpoiAs) - OM)) ■ Note that by 
symmetry G'poZ,iy(s', s) = GpoiA^',s). From Eqs. I11I9I 
we readily obtain for the full 4-block model 

Gpoi,x{s',s) ^ ^ / / (isids2cos((7o(s2 - Si)) 
^ Jo Jo 



exp 



lo Jo 

\S2 - Si 



(12) 



where for a free microtubule (not angularly constrained 
at the ends) we have to integrate over the arbitrary initial 
angle ip € [0, 27r] . The variance can now be computed and 
we obtain : 



2^ P^g f , 
Vpoiis -s) = ^^||s 



i-eq'o 



i+eq'o 

-2gfoCsin(go |s' - s|))} 



{{l-qie) cos(go(s'-s)) 

(13) 



with £, = £,b. From Eq. [13] the persistence length Ipoi 
and thus Ip can be deduced. The polymorphic variance 
Eq. [13] has the following generic asymptotic behaviors 
: At large distance s » ^, the variance Vpoi{s) ~ 
2k\P'^ i-i-l^g^ ^ scales linearly with s indicating that 9 poi 
is performing a simple (angular) random walk. At such a 
scale the microtubule looses its "coherent nature" and is 
replaced by a collection of uncorrelated segments. So not 
surprisingly we recover the classical results of a semiflex- 
ible chain again. In this asymptotic regime the effective 
persistence length reaches saturation with a renormalized 

constant value Ip (00) = 1/ + Z^^^ where 



Ipoi (00) 



l + eq'o 



(14) 



At short distance such that s <<min(^, g^""^) the vari- 
ance has a quadratic behavior Vpoi[s) « k\P^s^. In 
this regime, the polymorphic fluctuations are completely 
dominated by purely "classical" semiflexible chain fluc- 
tuations and lp{Q) — Ib- Starting from this value of Ip 
at s = 0, polymorphic fluctuations begin to contribute 
reducing the persistence length Ip (s) ~ {-f^ + "^^^ ^ )^"'"- 
This behavior is universal and independent of ^ and Qq^ ■ 

In the intermediate regime with s of the order of ^, 
the behavior of Vpoi{s) depends on the value of ^(70 - a- 
kind of "helix coherence" parameter. When S^qo < 1 (low 
helix coherence) oscillations in Eq. [13] are damped and 
the persistence length lp{s) is monotonously deceasing 
from Ib to the constant Ip {00) with Ipoi (00) « {£,i^iP^)~^ 
obtained from Eq. [14] (cf. Fig. 4) . This corresponds also 
to the situation of no internal twist go = 0. 

In the most interesting (high helical coherence) regime 
^qo > 1, the behavior of Vpoi{s) is a combination of two 
effects: an oscillation with wave length A = "^TrqQ^ orig- 
inating from the helical nature of the microtubule, but 
which is now damped by the presence of thermally in- 
duced defects (due to a finite £,) reducing the coherence 
of the helix and enhancing the linearly growing behavior 
(the random walk) . As a consequence for high helical co- 
herence Ip (s) displays three different regimes (cf Fig. 4): 

I. In the limit of very short , Ip (s) « (j^ + '^^^ '' )~"^ 
that attains a global minimum at Smin ~ ^Qo^ 

II. For intermediate length values Smin < s < ^, the 
total persistence length displays a non-monotonic oscilla- 
tory behavior of period A with damped amplitude around 
a nearly linearly growing average reflecting the polymor- 
phic fluctuation of the helix. 

III. For distances s >> ^, the oscillations in Eq. [T3] 
are completely damped as the helix forgets its "coherent 
nature" on these scales and reaches the limiting value 
i; (00) = 1/ (l-^i + l^^^ with Ipol (00) given by Eq. M 



HIGH COOPERATIVITY LIMIT AND UNIFORM 
STATES 

In the previous paragraphs we considered the ther- 
modynamic limit L — Nb — !■ 00 for which the correla- 
tion length ^ is always smaller that the total length L 
[ll|. It also interesting to consider the opposite regime 
^ >> L. In this case of large cooperativity we may 
take formally the limit J — > 00 . Consequently the 
two states variables ct* become uniform all along their 
respective block axes i — 1..4. The decoupling of the 
full MT lattice into 2 independent 2x2 sub-blocks 
implies the equality of the energy and entropy of the 
two independent 2-sub-blocks sets so that wc can fo- 
cus on only a single sub-lattice made of two blocks, 
say 1 and 2. The corresponding polymorphic energy is 
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FIG. 4: Representative plots of the angular persistence length 
Ip versus distance along contour s for intrinsic elastic bending 
stiffnes Ib = 10/im (chosen arbitrary and small for pictorial 
convenience) and different values of scaled correlation length 



E = [(§7 - 1) (ct^ + cr^) + 20-1(72] and its parti- 

tion function reads now Z = l + 2cxp ^—2(7 — B)L/bj + 

exp(— 47L/6). For two blocks there are then only 4 
possible states with different probabilities of realization 
p(CTi,cr2) = exp{-l3E)/Z that are: p(0,0) = 1/Z for 
the straight unstrained state, p(l, 1) = exp(— 47L/6) /Z 
for the straight prestrained state and finally p{0, 1) ~ 

p{l,0) = exp ^—2(7 — B)L/bj jZ for the two polymor- 
phic helical states. The average value of the spin and the 
correlation functions are given respectively by (f'^) = 
/9(1,0) + for k = 1..4 and (fiV^) = (a^a^) = 

Going back to the full lattice model made of four blocks 
the polymorphic order parameter is then given P^{L) — 
2{{a^)-{a^a^) + {a^)-{a^a*)) =M1,0). The angular 
variance at position s and s' can be easily computed and 



we find V{s' — s) 



(1 — COS (go(s' ~ s))) which 



obviously corresponds to the limit of infinite ^ of Eq. [T3] 
but with the additional factor P^{L) responsible for the 
length dependence of the variance. The weight of the 
different possible configurations depends on 7 and the 
discrimination between them is more pronounced with 
growing length. This can be easily understood from the 
polymorphic entropy Eq. which remarkably is non ex- 
tensive. Note that E^r is always twice the entropy of 
the individual 2-block-sub-lattices. For the asymptotic 
limiting case of short L — all 4 x 4 configurations of 
the full lattice have the same probability and therefore 
Eo-(O) = 2^5 In 4. For larger L, the entropy will continu- 
ously change to reach a constant value depending on the 
dominant configurations that are selected by 7. For val- 
ues of I7I > B, we see on Fig. 5 that P{L) goes to zero 
with growing values of L. In this regime there is only one 
configuration (the same for the two sublattices) p(0, 0) 



or p(l, 1) depending on the sign of 7 and Ecr(i) — ^ 0. 
Interestingly for 7 = ±S, there is an additional config- 
uration for a given sub-lattice with the same probability 
p(0, 1) = p(l, 1) = p(0, 0) = 1/Z, so that for the 4-blocks 
we have 3x3 configurations and ^^(L) — > 2kB In 3. This 
coexistence leads to a state which is the superposition 
between straight and curved states with different prob- 
ability such that P^iL) 



that is 



3+cxp(-4BL/6) ~ 3 

(quasi) length independent. When I7I is very close to B 
, the evolution of Eo-(i) versus L, can be slow and there- 
fore there is a regime of lengths where the the entropy is 
non-extensive and where the straight states have compa- 
rable free energy with the helical one. In this regime the 
microtubule will fluctuate between its almost degenerate 
straight and (two) curved states. 

For I7I < B, there are for a given sub-lattice two dom- 
inant degenerate configurations (0,1) and (1,0), so that 
the number of configuration for the full lattice model is 
2x2 and Eo-(i) — 2fcBln2 with growing L. In this 
regime the mean curvature is built up progressively with 
the length L, to reach -faster for smaller value of I7I - its 
maximum value P = a/2- Here also when I7I is very close 
to B the realization of the curved state can be slow with 
growing L. 




FIG. 5: Polymorphic entropies E (four-block model) versus 
the microtubule's length L for different values of conforma- 
tional bias 7 and stiffness B — 1. 



ANGULAR POWER SPECTRUM 

Other information on the equilibrium properties of 

the microtubule can be deduced from the analysis of 

its Fourier mode distribution. To do so we express the 

microtubule shape 9 (s) as a superposition of Fourier 

modes 6 (s) = 
,+00 



wave vector q„ 



'(0) + v^E„ri^nCos(g„s) + 

. sin(g„s) -|- {9j;{L) - 9x{0j) s/L with the 
and where we add the "ramp term" 



2n7r 
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P {a^) 




FIG. 6: Plots of the polymorphic polarization P (four-block 
model) versus the microtubule's length L for different values 
of 7 and with B — 1. We observe that microtubule curvature 
can be build up progressively with growing MTs. 

{9x{L) — 9x{0)) 3 / L, as 9 (s) is not a periodic function 

in general, i.e., 9 (L) ^ 9 (0) [l^. Decomposing the 

shape as "^(s) — 0poi{s) + 9 ei{s) we readily see that 
the equilibrium square of the amplitude of the elastic 

Fourier modes Geiiqn) = {'^li,n) = (^li.n) = 
has the usual thermally driven elastic bending spectrum 
of a wormlike chain. On the other hand from Eq. [12] 
the polymorphic mode distribution Gpoi{q) in the limit 
L large (with qn ^ q continuous) can be deduced : 

(15) 

As shown in Fig. 6, when the correlation length ^ is of 
order of the helix pitch A = 27r/go, the function Gpoi{q) 
departs very much from the wormlikc chain bending spec- 
trum Ge/ {q) and displays a non- monotonous shape with a 
peak at go- This peak is more pronounced and sharper for 
larger value of (helix more coherent) and disappears 
when 1- This is the main Fourier mode signature of 

the -polymorphic fluctuating- helicoidal nature of the mi- 
crotubule. For long wavelength modes q « qo, the spec- 
trum has a semiflexible like behavior Gpoi [p) ~ j — ^^oo)q'^ 
reflecting the incoherence of the helix at this scale, with 
the persistence length given by Eq. [TH For short wave- 
length modes q » qo, the behavior becomes quarticaly 
decreasing Gpoi{p) ~ and polymorphic fluctuations 
are strongly damped at short distance in agreement with 
high cooperativity at this scale. 

Quantitative comparison with experiments at this 
point seems difficult, because Fourier analysis of micro- 
tubules in-vitro was mainly performed on microtubules 
strongly confined into an almost 2-D geometry. This ex- 
perimental necessity (originating from optical microscopy 



FIG. 7: Log-Log plot of the polymorphic Fourier mode dis- 
tribution Gpoi{q) for different values of (,qo- Besides a charac- 
teristic peak at qo one observes a cross-over between a small 
q regime with scaling q^^ and a large q regime with a 
scaling. 

sample flatness requirements) is expected to induce arti- 
facts as in this case, within our helical model, the micro- 
tubules would likely not fully equilibrate between their 
different equivalent curved states as large barriers will 
appear under confinement. In agreement with this ar- 
gument the experimental Fourier modes of microtubules 
systematically reveal some type of " frozen- in" curvature 
which is much larger than the fluctuations around it [isj 
- in particular for longer wavelength modes. This is in 
so far striking as no clear argument for large system- 
atic built-in curvature in microtubules (other than the 
here described) is evident, indicating that this curvature 
is in fact the very character of the microtubule elastic- 
ity. Similarly in vivo, it was found that growing micro- 
tubules have large scale frozen- in curvature [14]. This 
phenomenon was interpreted in terms of a random inter- 
action of the growing microtubule tip with the surround- 
ing cytoskeleton. Instead it could be that the curvature 
has an intrinsic origin (like in the present model) but is 
not allowed to fully equilibrate due to constraints present 
in the cytoplasm. 



CONCLUSION AND OUTLOOK 

We have investigated the thermodynamics of a toy 
model for cooperatively switching (polymorphic) bioflla- 
ments, whose paradigm example is believed to be the mi- 
crotubule. The polymorphism of the microtubule's sub- 
units and the short range cooperative interaction leads 
to properties that are very different from usual biofila- 
ments. In particular the ground state itself is not a single 
fixed shape but is found to consist of a set of degener- 
ate helicoidal configurations. At finite temperature, the 
fluctuations of the subunits create defects that smoothen 
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the perfect helix and can even destroy it on very long 
scale where the filament retrieves a typical worm-like 
chain behavior. Another peculiar characteristic of the 
polymorphic tube is the length dependence of the persis- 
tence length. The Fourier spectrum turns out also to be 
markedly different from classical biofilaments. While it 
has a typical semiflexible chain behavior at long scale the 
Fourier modes amplitude are enhanced for modes around 
the characteristic scale of the helix. At very short scales 
the modes are strongly dampened due to the strong co- 
operative dimer interaction. We believe that the Fourier 
spectrum will be of some value to discriminate between 
different models in the analysis of future experiments (as 
in [J]) where confinement is absent and microtubules can 
thus rearrange and equilibrate. 

The model developed here was inspired and tailored 
to microtubules. However it can be adapted to other 
polymorphic filaments with similar cross-sectional sym- 
metry (circular rod or tube). Polymorphic sheets, rods 
and tubes are not uncommon in nature with^ the best 
documented example the bacterial flagellum |l5j. In the 
latter case the polymorphic states are likely very rigid 
(or coherent in our terminology) on thermal energy scales 
with correlation lengths ^ larger than filament's length L. 
Rearrangements between equivalent lattice states could 
then face high barriers in the fiagellum and suppress the 
here described polymorphic fiuctuations on measurement 
time scales. However other filaments might be even bet- 
ter candidates for soft polymorphic tubes (or rods) and 
deserve a closer investigation. Observations of unusual 
fiaring-up of peculiar helical modes in bacterial pili (cf. 
Fig. 1 in [161] ) as well as in actin filaments cooperatively 
interacting with drugs like cofilin (cf. Fig. Ic in (l7| ) 
visually suggest the presence of soft polymorphic mode 
dynamics. 
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